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2. 


0-  Introduction .  Consider  s  (>  2)  linear  regression 


models 


(0.1)  Xki  ak  +  6kcki  +  Zki  .  i  1 . nk  ;  k  1 . S 

where,  for  each  k  =  l,...,s,  a,  is  the  (unknown)  intercept, 
(0.2)  8k  =  (8kl, • . . , 8kq) 

is  a  q-dimensional  vector  of  unknown  regression  parameters, 


(0,3)  Ski  (ckli' - Ckqi} 

is  a  q-dimensional  vector  of  known  regression  constants  for  each 
i  =  /  and  the  Z^  are  all  independent  (error)  random 

variables  with  the  same  (but  unknown)  continuous  distribution 

(0*4)  F (x)  —  P( Z^ ^  <x)  ,  k  -  1,  . .  .  , s  7  i  =  1 ,  •  •  •  , n^  • 

A  problem  of  interest  is  that  of  testing  whether  the  s  regres¬ 


sion  surfaces  are  parallel  to  one  another,  i.e.. 


(0.5) 


H  :  ft.  =  .  .  .  =  6  =  S (unknown) 

0-1  -  s 


Accession  lor 

NTIS  CM*! 

Hire  ?;.»<  ;  ; 

ire.’.,:*  . 


H  :  B,  /  6  •  for  some  1  <  k  ^  j  '  s 

—  K.  —  1 


For  the  special  case  q  =  1  ,  i.e.,  testing  the  parallelism  of  sev¬ 
eral  regression  lines.  Sen  (1969)  has  proposed  a  class  of  rank  order 
tests.  In  the  present  paper  we  study  the  problem  in  the  general 
case  q  >  1  .  Preliminary  notations  and  assumptions  are  given  in 


Alii  i'vS  ,!h  t_ 
N  *7  I  •  - 


!.■  M  .’NTIFIC  KKShAKOU  (AFSCJ 


“  ''  and  Is 

iaa  a: h  i^u-12  (?b). 


A.  i). 

Tt'o’hnicai  iMVinritlon  OTficer 


section  1.  In  section  2  a  class  of  asymptotically  dist r ibution- 
iree  aligned  rank-oroer  tests  are  proposed.  These  are  univariate 
counterparts  of  a  nu \ t Lvariat e  problem  briefly  mentioned  in  Sen 
and  Puri  (1977)  but  not  solved.  The  asymptotic  distribution  of 
t he  test  statistics  is  derived  in  section  3.  In  section  4  we 
derive  the  asymptotic  relative  efficiency  of  the  proposed  tests 
with  respect  to  the  general  livelihood  ratio  test  of  the  same  pro¬ 
blem.  Finally,  in  section  5  asymptotic  optimality  in  the  sense 
o'  w .  ■  ’  d  {  l  9  4  1  is  d  i  s  o  ur  '*  ed  . 


!■  •  Preliminary  Notations  and  Assumptions.  For  each 
k  =  l ,  .  .  .  ,  s  let. 


(1.1) 


J,  1  n. . 5kqn  / 


where 


(1.2) 


-  _  -i 

Ckmn^  r‘k  3^  ckmi  '  111  *  1 . 01  * 


We  assume  that  the  cr  *  q  symmetric  matrices 


(1.3) 


M 


kn, 


nk 

V 

i-1 


[  c 


:ki 


c  ) 
kn, 
k 


/  k  1  f  •  i  s 


are  positive  definite  and  that  the  limiting  matrices 


(1.4) 


-  lim  n"1  Mkn  ,  k  =  1 . s 

k 


exist  and  are  positive  definite. 


Simplifying  some  of  Jureckova's 


4. 


(1971)  conditions  on  the  regression  constants,  we  also  assume 
that  each  c^  can  be  expressed  as  a  difference 


(1.5) 


Ski  Ski  (1)  cki (2)  ' 

cki(j)  =  (ckli ( j ) . Ckqi(j)} 


where,  for  each  X  =  ,  m  =  1 and  j  ~  1,2,  c 

Is  nondecreasing  in  i  ,  and  the  c^mi(j)S  satisfy 


kmi  ( j) 


(1.6) 


-1  _  2 
lira  n.  max  c.  .  ,  .  v  -  c,  ,  .  N  I  -  0  , 
nk..„  X  Xmi(j)  to^U) 


-  -If 

'kmn^(j)  nk  ckmi(j) 


(1.7)  lim  n,  1  T  'c,  ,  -  c,  .  e  (0,») 

n  k  i4i  ^i(j)  kmnk(j)J 

X 


which  together  imply  the  Noether  condition 


(1*8) 

lirn  f  max 
nk-  l<i  nk 


kmi ( j ) 


2  /  K  _  2 

kmnk(j)  '  [Ckmi(j)  "  ckmnk(j)1  *  =  0 


We  denote  the  total  (combined)  sample  size  by 


(1.9) 


»  =  1 


and  assume  that  the  limits 


5 


(1.10) 

rk  ~ 

lim  (n  /N)  ,  k  =  1,  . 

exist  and 

satisfy 

(1.11) 

r  < 
o 

rk  5  1  -  ro  •  *  *  + 

for  some 

0  <  r  < 
o 

1/s  .  Thus  we  have 

(1.12) 

lim 

V-vco 

=  °°,  k  =  1 ,  . . . ,  s 

and 

(1.13) 

S 

l 

"k  =  1  • 

k=l 

and  the  matrices 


(1.14)  M*  »  lim  N~^K.  =  lim  (n./N)M.  =  r.M.  ,  k  =  1 _ ,s 

K  Knx  Nh-oo  k  k  k  k 

are  symmetric  and  positive  definite. 

For  each  positive  integer  n  ,  let  the  scores  a  (l),...,a  (n) 

n  n 

be  generated  by  a  non-constant  and  square  integrable  function  \fj 
on  (0,1)  according  to  one  of  the  following  two  ways: 

(1.15)  an(i)  =  tj>l!i/(n  +  1)1,  i  =  l,...,n 

or 

(1.16)  a  (i)  =  Efti(U  .)!  ,  i  =  l,...,n 

n  ni 

where  U  .  £  ...  ■?  U  _  are  the  order  statistics  of  a  random  sample 
ru  nn  ^ 

of  si2e  n  from  the  uniform  distribution  over  (0,1).  We  assume 
that  4>  can  be  expressed  as  the  difference  t|»  s  i|;.  ^  of  two 


6. 


non-decreasing  and  absolutely  continuous  functions  and 

on  (0,1)  .  Let 


(1.17) 


X(ifr) 


[4>(u)  -  $]2  du}*5 


4' 


4-'  (u)  du  . 


Thus  we  have  0  <  A  (4‘)  <  00  . 

We  assume  that  the  underlying  distribution  function  F  has 
an  absolutely  continuous  density  f  =  F*  with  finite  positive 
Fisher  information 


(1.18) 


0  <  1(f) 


[f ’ (x)/f (x) 1 2dF (x) 


CO 


We  note  that 


(1.19)  1(f) 
where 

(1.20)  <frf(u) 


•f 


[  X  (4>f)  !2  =  J  C<>f  (u)  l2du 


-f'[F"1(u)]/f[F“1(u)  ]  ,  u  e  (0,1) 


with 

(1.21) 


4>f(u)du  = 


0 


2.  The  Proposed  Rank-order  Tests .  For  b 
and  k  *  1, . . . ,  s  ,  let 


(b,  - - b  )  c  lRq 

i  q 


(2.1) 


Rkink<?)  '  the  ranX  °f  Xki 
Xkl  “  bckl . Xknk  "  b?knk 


-  be,  .  among 
»  *  Kl 

in  the  ascending  order  of 


magnitude 


7 


nk 

12'21  Skmr.  (?>  '  .1.  (ctoni  "  5kmn>  anvl  Rkin„ <?' 5 '  m  =  1 . » 

k  i~l  k  k  k 


where  a  (l),...,a  (n.)  are  generated  according  to  (1.15) 

nk  nk  K 

or  (1.16)  (with  n  replaced  by  nv) , 


(2'3)  Sknk(b)  ~  (S*lnk(b) - ,Skqnk(b))  ' 


and  define 


S 


(2.4)  Ss(b)  .  ^  Stojt(b)  -  (SN1(W . SmM)  . 

Let 


(2.5) 


B(N)  =  {*>  e  *  ?  I  S  (b)  >  =  minimum} 

v "  '  ~  m=l  l'‘  ‘ 


and  choose  one  element 


(2.6) 


~N  €  B(N) 


as  an  estimate  of  8  .  Define  the  s  vectors  of  aligned  rank 
statistics 


(2.7) 

and  let 

(2.8) 

where 

(2.9) 


f  Nk  "  ?kn^  '  k  1 . s 


»„  "  tx'1  l  l  <i>  -  *n  1“) 

LN  k~l  i  =  l  nk  nk 


s  k 


n. 


-1 


a  -  n.  T  a„  (i)  ,  k  =  l,....s 

nk  k  iii  "k 


8. 


Then  a  class  of  aligned  rank-order  tests  of  (0.5),  each  deter¬ 
mined  by  a  score -generating  function  ,  can  be  based  on  the 
statistics 


A  + 


(2.10)  Q[)  =  I'2  Sm 


whose  asymptotic  distribution  under  Hq  is  given  by  Theorem 
2.1,  which  in  turn  follows  from  Theorem  3.1  (see  section  3). 


Theorem  2.1.  Under  Hq  ,  QN  has  asymptotically  the  (central) 

2 

chi-square  distribution  *(s_i)q  with  (s  -  l)q  degrees  of  free¬ 
dom. 

For  0  <  e  <  1,  let  X^  . .  be  the  upper  100e*  point 

V s  —  X /  C[  /  & 

2 

of  the  X.  . .  distribution.  Then  for  large  N  we  have  the 
(s-l)q 

following  asymptotically  distribution-free  test  of  approximately 
size  e  : 

(2.11)  Reject  Hq  (in  favor  of  H  )  if  and  only  if 
°N  S  X(s-l)q,e 


3.  Asymptotic  Distribution  of  the  Test  Statistics.  Consider 
the  sequence  of  hypotheses 


(3.1)  hn  :  6k  3kN  6  +  N  **bk  ,  k  1, 


where  the  s  vectors  b*  e  ]R^  ,  k  =  l,...,s  are  such  that 


l  £k  *k  =  9 

k_l  ~K  K 


(3.2) 


9. 


10 


results  (see  Theorems  3.1  and  4.1,  and  Lemmas  4.1  and  4.5) 

A 

the  distribution  (denoted  by  V  )  of  n(?(3  -  6  )  is 

K  ^  kn^  ~  K 

asymptotically  normal  (denoted  by  N^) ,  i.e., 

(3.8)  Prn^(Skn  -  6^)]-^  V  (0, [ A ty) /y (ip , f ) )2m“1)  (k  =  1, 
k 

and 


(3.9) 

Similarly  we  have 

(3.10) 


n  %  (6,  )  =  o  (1)  (k  =  1 . s) 

k  - kn^  ~  knk  p 


'’"’VlV  -  °p<1’ 


We  need  the  following  lemmas  to  prove  Theorem  3.1. 


Lemma  3.2.  For  each  k  =  l,...,s  we  have 


-hi 


(3-U>  "'Idk  -  ’'('i”£,N!5(?kr,k  -  5d1Mk  +  V1' 


Proof .  By  Theorem  3.1  of  Jureckova  (1971),  for  each 
k  =  1 , . . . , s  we  have 


(3.12) 


and 


2Sv"  )  nk*?knk(§k) 


‘k  -knk'lknk 


"  ^'f)nk(Bknv  ~  Bk)Mk  +  °p(1) 

K 


nkH?kn)((?N)  =  "k^Hn^lhc* 


r(#.f)n^(SN  -  6k)Mk  +  op(l)  . 


(3.13) 


11. 


Subtracting  (3.12)  from  (3.13)  and  using  (2.7)  and  (3.9),  we 
have 


(3.14)  nk^SNk  Y  (v' .  f ) kn^  ^N^Mk  + 

l 

Multiplying  both  sides  of (3. 14)  by  (n^/N)  ? 
we  obtain  (3.11) . 


°pC)  . 

and  using 


(1.14) , 


For  later  use  we  also  define  the  q  v  q  mati  ix 


(3.15) 


D  = 


s 

r1 

l 

k=l 


M 


*  __ 


k 


=  (d  c)  „ 

m.t  m,  x= I 


•  a 


which,  being  a  sura  of  symmetric  and  positive  definite  matrices, 
is  itself  symmetric  and  positive  definite  and  hence  has  a  symmetric 
inverse 


(3.16) 


A  =  D 


-1 


^a£m^  £,m=l. 


<  g 


Thus  we  have 


(3.17) 


s  s 

r  . 


DA  =  AD  =  l  ML  A  =  A  [  Mv 
k=l  K  k=l 


where  I  is  the  q  x  q  identity  matrix. 

q 

Notation.  Let  {U  }  and  {Vn>  be  two  sequences  of  random 
vectors  of  the  same  dimension.  Then 


(3.18)  un  "  Vn  if  and  only  if  un  -  =  o  (1)  . 


Lemma  3.3. 


N^S 


N 


N 


S  A 

l  6 

k=l 


knk^ 


(3.19) 


13. 


where 


14 


(3.27) 


(t 


kl' 


■  V  < *q 


k  1  p  •  •  •  t  S 


and  let 


(3.28) 

Then  oy  (3.19)  we  have 


u  =  l  t.  . 

k=l 


'l  ^  A  A 

T  t '  -  7  N  *  (  9  -  '•?  )  t/ 

N  k^i  -knk  'N  k 


(3.29) 


T  tr  -  N*73  u‘ 

k-1  ~krVk  N 


S  1,A  1.  S  A 

I  N2Skn  fck  -  N*  l  &kn  MkAu' 
k=l  ~knk'K  k=l  knk  k  ' 


A  A 

By  making  the  substitution  s,  =  ( 8,  -  P)  +  8  on  the  right- 

"knk  ~Knk 

hand  side  of  ~  in  (3.29)  and  then  making  cancellation  (using 
(3.17))/  we  have 


(3.30) 


-N 


■ 31  (t*  ■  ^flu'’ 


Now  by  (3.21)  and  the  symmetry  of  and  A  ,  under  the 

random  variable 


N  (6kn  ~  V  (5k  “  MkA^} 
k 


is  asymptot ically  normal  with  mean  b*(t^  -  M* Air  )  and  variance 
l\(i(-)/y(i,f)  !2(tk  -  uAM*)M*-1(t;  -  M*Au')  . 


V  k 


So,  by  independence  of  the  fikn  '  s  ,  under  the  right-hand 


15. 


side  of  (3.30)  is  asymptotically  normal  with  variance 

[  \  (.  )  /  ■  ((  ,  f )  '  2c2  where 


s 

(3.31)  c2  =  l  (tk  -  uAM^)M*-1(tk  -  M*Au') 

k=l 

and  mean 


I. 

k=l 


s 

l 

k=l 


b*t: 
-k~  k 


(  l  b*M*)Au' 
k~l  ' k  k  ' 


b*t" 


the  last  equality  in  (3.32)  being  a  consequence  of  (3.2). 
Expanding  the  right-hand  side  of  (3.31)  and  using  (3.17),  we  have 


=2  =  Jj  hCY*  -  Jj  hAu' 


=  h  J'Vj'1  -  A,h 


=  tJt" 


qr*  v 

Thus,  for  any  t  e  13"  ,  TNt  under  has  asymptotically  a 

(possibly  degenerate)  normal  distribution  with  mean  b*t'  and 
variance  [  \  (ij>)  /y  ,  f )  ]2tJt'  .  It  follows  that 

P(Tn|Hn)  H  lb*,  [A(^)/y(iJ>,f)  ]2J)  . 


(3.33) 


16. 


Lemma  3._6.  Under  H..  ,  the  random  variable 

N  ‘  - - 


(3.34)  q*  -  J  <5to  -  -  V' 


is  asymptotically  ^  (S-l)  q  * 


Proof.  Let 


(3.35)  Yn  =  rY(^f)/X^)]TN  =  (Yn1 . YNs) 


where 


(3*36)  Ljk  =  [Y 


(^.f)/x(^) ]N^(ekn  -  tN)  ,  k  =  i . s 

k 


Then  by  (3.33)  we  have 


(3.37) 


P(vn1hn)  -*■  ysq([r(U)A(t)lb  ,j> 


Define  the  (sa)  *  (sq)  symmetric  matrix 


(3.38)  K  {Kj<j)x,  j=i,  . .  .  ,s  (6kjMk)k,j  =  l . s  ' 

Then  (3.34)  and  (3.3)  can  be  rewritten  respectively  as 


(3.39) 


°N  L  YNk^^Nk  YNKYN 
k=l 


(3.40) 


Ag  =  {  r  Y  (4»,  f )  /X  (i|j)  ]b  }K{[7  (t|),  f)/X  (\|))  ]b*}'  . 


So,  to  prove  Lemma  3.6,  it  suffices  to  show  that  KJ,  or  equiv¬ 
alently  its  transpose 


(3.41) 


W  =  JK  , 


17. 


is  id* mpotert  with  tr<ice  equal  to  (s  -  l)q  and  that 


3.42) 


b*KJKb*'  -  b*Kb* ' 


(see  (!°7l)  ,  Corollary  2s.  1)  .  By  direct  computation 

we  have 


(3.43^ 


W  *  fW.  .K  .  , 

kj  k, : -1 ,  .  .  .  ,  s 


(3.44: 


W.  .  -  ; .  .  I  -  AM  ,  k,  j 

kj  ki  q  J 


3y  further  computation  and  (3.17)  we  have  -  W  .  On  the  other 


KJK  =  KW  =  (5  Mj  -  MkAM))k 


and  so 


b*KJKb* 


(3.45) 


V  b*M£b* ' 
k  vvA 

k~l 


=  b*Kb 


(  "  b*M*'A  •  M*b* 
v  - 1  K  k  j  - 1  J  J 


where  the  last  equality  in  (3.45)  follows  from  (*.2) 


It  remains  to  compute  the  trace  of  W  .  Let. 


«3-46>  Mk  "  (ckn,«>».  (=1 . q  ’ 


k  =  1 


Then  by  (3.15)  we  have 


(3.47) 


dm<L  “ 


,  m,  l  =  1  , 


*nd  v  [  * .  ’  *■' )  and  (  * .  L  7 }  we  have 


1 


1  z  •  .  .  /  cj 


*ow  :  .'T  k  *  l 


by  (3.44)  we  have 


W  -  I  -  AM, 
V  *■  <  t  v 


a  ,  e.  .  .  )  / 

•  r  K  Tt  •  -  ,  »  '  1 


f  w .  •  ;  *  •  ?  t 


.pr  km  * 


l  S  .  4  ■.* ' 


t  I  <Vi) 


'■  t  'wkk'  r  -  q  ~  l ) q 


'hus  !  .ep'irut  3  .  *>  v  ^  estab^ 


'r°c>f  of  rheoren  3 .  i  .  Hy  herrma  3,6  it  suffices  to  show  that 


M .  •  ' 


by  (l.UM,  (l  \r>)  -  11.17'  and  (2.0)  we  have 


1  *  *r>  r>  ^ 


and  so  by  (<..8*.  (1.10)  and  '1.13)  we  have 


(301) 


1  un  \  -  X  U) 

N  »«’  N 


It  follows  from  (3.11)  and  (3.36)  that 


(3.52) 


,NXNr'‘?'Nk  ’bX  ■  k  *  1 


19. 


Now  by  (1.14)  we  have 


(3.53)  lim  NMr^  =  M*_1  ,  k  =  1 . s 

N-vcw  nk 

It  follows  from  the  symmetry  of  that 


(3.54) 


-2A  -1  A. 

XN  SNkMki  SNk 
k 


-  y  m*y" 

Nk  kCNk 


Summing  up  both  sides  of  (3.54)  over  k  =  and  using 

(2.10)  ard  (3.39) ,  we  obtain  (3.50).  Thus  Theorem  3.1  is  proved. 


4.  Asymptotic  Efficiency.  Using  (3.~M,  we  rewrite  (0.1)  as 


(4.1) 

where 

(4.2) 

(4.3) 

and 


X,  =  +  e,  Ci  +  Z, 

~kn^  k~  -k  k  ^kn^ 


k  =-  ]  ,  .  .  .  ,  s 


n. 


1  -  (1 . 1)  €  * 


n 


Z.  “  ( Z.  .  /  •  •  •  t  Z.  ) 

kn^  kl  kn^ 


ck  *  (ckl . 5knk> 


(4.4) 

is  a  q  v  matrix. 

Let 


(4.5) 

(4.6) 


X  =  ( X .  . X  ) 

- N  V  in,  -  sn 

1  s 


- N  ~  in,  - sn_ 

1  s 


20. 


(4.7)  a  =  (ctw...,a)  , 
-is 

(4.8)  e  =  (a, . . .  ,e  )  . 
and  let 


(4.9)  E,  =  (0: - r  _ ,0')'  ,  k  -  1 . s 

k  -  ~nk 

be  the  s  x  n.  matrix  with  1  as  the  k-th  row  and  all  the 
K  ~nk 

other  rows  being  0  -  Then  the  s  linear  models  in  (4.1)  can 
be  combined  into  one  linear  model 


(4.10) 


X  S!  PC  +  Z 

Cn  :tn 


where 


(4.11) 


*  _ 


E1  E2 
c1  0 

0  c„ 


0 

0 


.  0  0  .  .  .  c 

*■  s 


is  an  r s (q  +  1)  j  x  N  matrix.  The  parameter  space  for  0  is 
the  s(q  +  1) -dimensional  Euclidean  space 


(4.12) 


and  Hq  can  be  expressed  as 


(4.13)  H  :  0  t  GL  =  {(a,b1(...,b  )  e  ft  :  b,  =  . .  .  =  b  }  . 

o  o  -  - 1  - s  -1  ~s 


The  likelihood  ratio  test  of  (0.5)  rejects  Hq  (in  favor  of  H  ) 


J 


21. 


if  the  likelihood  ratio 


s  nk 


‘lN  SUP-1  Jl,  f{\±  -  ak  -  b5ki>  S  ak  €  R  #  k  =  1 . s;  b  c  ]Rq}/ 


k'l  l-l 


s  nk 


k-  1  i-.l 


f(xki  -  ak  -  bk?ki'  =  «k  «  *  .  bk  «  *<»  .  k  -  1 . s! 


is  small,  or  equivalently  if 

(4.14)  L =  -2  loq  A„ 

v  N 

is  large.  Her*'  f  (or  equivalently  F  )  is  assumed  to  be  known. 
Under  Assumptions  I  -  V  and  VII  of  Wald  (1943),  but  no  assumption 
concerning  the  shape  of  F  ,  the  asymptotic  distribution  of  L 

N 

under  is  qiven  by  Theorem  4.1,  which  will  be  proved  later  in 

this  section. 

I  he  or  err  4_.  _1.  Under  H  ,  L  is  asymptotically  X  2  (AT  )  with 

^  l  s  —  1 ;  q  L  — 

al  •  iwX  bXbk"  • 

To  compare  the  proposed  rank-order  tests  with  the  likelihood 
ratio  test,  we  make  the  additional  assumption  that 

(4.16)  b*  /  0  for  some  1  <  k  *  s  , 

which  makes  the  right-hand  side  of  (4.15)  strictly  positive.  Com- 
bimm  Theorems  3.1  and  4.1,  we  have  the  asymptotic  relative  effi¬ 


ciency  . 


22. 


23. 


Let 


(4.20' 


r  ( i )  <§o'h)  ^l"* 

■■Zj  =  (z 


^s+q^ 


^  (2)  "  (§2 


s+q+1 '  *  * ‘ ' -s (q+1) ) 


Then  can  be  expressed  as 


(4.21) 


§( 2)  (e)  h  o  . 


Clearly  £  is  a  homeomorphism  and,  with  the  identification 


(4.22)  e  =  (a.?2 - -  8S)  =  ( 0 ,  .  .. ,  0g  (  +1)  ) 

has  a  positive  Jacobian  det Og/98)  not  depending  on  0  ;  more¬ 
over,  the  first  two  partial  derivatives  of  E-.  (0) . E  ,  .  (6) 

1  -  Ss  (q+1)  ■- 

are  uniformly  continuous  and  bounded  functions  of  e  .  Indeed, 
the  inverse  @  =  E~^  of  £  is  given  by 


(4.23) 


e(E)  =  (§o'fl'§l*S2 . h+^s] 


with  Jacobian  matrix 


(4.24) 


M  =  (30/3E)  = 


0 

I 


0  , 

0  . 
q 

cq  Iq‘ 


o 

o 


o  i  o  ...  i 
q  q 


(whose  determinant  is  equal  to  unity) . 

Now  consider  the  ’s(q  +  1)1  x  [s(q  +  1) ]  matrix 


24. 


(4.25) 


A  -  C*C* 
N  N  N 


F1  F2  •••  F, 


rl  U1 
0 


0 

D. 


.  0 

.  0 


F  0  0  .  .  .  D 


where 


(4.26) 


E  = 


^  5kjnk^  k, j  =  l , .  .  . ,  s 


(4.27) 


nk 

Dk  =  CkCk  =  i|1  ?ki ~ki  '  k  "  1 . s 


and 


(4.28)  FR  =  (O' _ 'nkSkn  . O')  -  k  =  1 . s 

k 

is  a  q  x  s  matrix  with  n]cc)cn  as  k-th  column  and  0 

K 

elsewhere.  By  routine  computation  we  have  the 
[s(q  +  1)]  *  '_s(q  +  1)1  matrix 


(4.29) 


AN  =  M  ANM  = 


^11 


N12 


N21  N22 


where 


(4.30) 


wll 


s 

l 

k=l  f 


s  s 

l  Fk  l  Dk 

k=l  K  k=l  K 


is  (s  +  q)  x  (s  +  q)  , 


25. 


(4.31) 

(4.32) 


and 


^N21  AN12 


(4.33) 


N22 


---  ( 


5,  .D.  ) 
kj  k 


k,j=2. 


We  note  that  by  assumption  is  positive  definite,  hence  so 

is  -  Consider  the  matrix 


(4.34) 


s 

V  y 

k-i  '  knk 


which  is  symmetric  and  positive  definite  and  hence  has  a  symmetric 
inverse 


(4.35)  G[J  - 

and  define  the  a  *  s  matrix 


(4.36)  C  =  (c  . c  ) 

N  -In,  - sn 

1  s 

Then  by  routine  computation  and  the  obvious  identity 


(4.371 


M. 


kn-. 


D.  -  n.  c.  c /  ,  k  =  1,  .  .  .  ,  s 

k  k-kn^-kn^ 


it  can  be  checked  that 


(4.38) 


\*-1 

mi 


cmgmcm+e 

N  N  N 


gncn 


-1 


■CNGN 


N 


26. 


with 


•  •  / 


(I)  -  (b* . b*)Aj(b* . b*) 


(4.46) 


T  b*M, 


-  '  l.  Sj'W'V  *  "kn 


ID” 

K=2  ~irkrVk  k=2  “  ~“k  "  k=2 


(II)  -  (b~ '  •  •  •  '  bg)  ^~1 '  *  *  '  '  ~  1^ 


(4.47) 


■2' 

s 


*  ^"kn^Vln^r 


(4.48) 


(III)  =  (II) 


and 


*  *  T~*  / ,  * 


(4.49) 


( IV)  —  (b^ ,  . . .  »bj_)  '  *  *  * '  ~1^ 

-  -  5l«lniVln^r  • 


By  (1.14)  and  (3.2)  we  have 


(4.50) 


And  by 


(4.51) 


“»  »'  J,  Sk"knv  = 

N-*»  k=2  k 


(3.15)  -  (3.16)  and  (4.34)  -  (4.35)  we  have 


lim  NG..  =  A 


N+co 


N 


It  follows  that 


lim  A*  =  T  b*M*b*'  and  so  lim  A 
m—  N  k=l  'k  lN 


=  A. 


The  proof  is  complete. 


28. 


5 .  Asymptotic 


-ity.  Let  and  be  non- 

-  NJ  N 


singular  square  matrices  of  orders  s  +  q  and  (s  -  l)q  / 
respectively ,  satisfying 


(5.1) 


■N1  =  Nil 


(5.2) 


bnbn  =  \ 


and  define  the  (s  +  q)  x  :  (s  -  l)q]  matrix 


(5.3) 


-  X-l, 


=  n!2 


Then  the  square  matrix 


(5.4) 


r  .  r 

N1  '  N2 


0  B. 


of  order  s(q  +  1)  is  nonsingular  and  satisfies 


(5.5) 


VT1^  *  ^(q  +  1) 


For  (o  =  (a,b,  ...,b)  c  SIq  and  c  >  0  define  the  surface 
(5.6)  S  (u> ,  c) 

-  (0  «  v  :  I(f)f;(2)  (6)Aj§(2)  (0)  "  =  c,  §(0)^'=  (a.b)!^} 


where 


(5.7) 


r  =  ( r  r  ) 
N  V  Nl'  N2; 


is  (s  +  q)  x  ! s(q  +  1) 1  .  Consider  the  transformation  of  ft 


29. 


(5.8)  r  -  ('w^ - -  -;s)  ~  . es)  =  [  Kf)  i^iejKj 

where 


(5.9) 


(«* ' !?* )  =  '  Kf)  J^5(0)r' 


and 


(5.10) 


(?2 . f*s>  =  '  1<f)  '  **  5  <  2 )  <!)Bn 


which  maps  ;j(u,c)  into 


(5.11) 


S  ( u  ,  c ) 


=  {  (/,  3*, 


★ 

?.>  < 


,  *  .  ★  . 
(a,  8^) 


K£)’'i<?-?>rsr  J2  ?!&'  "  c) 


For  Qq  <  P  and  p  >  0  let 


(5.12) 

P  (  0  Q  .  P  ) 

=  {8  r  Q:0,O^  c  S(w,c)  for  some  u  e  Pq  and  c>0,  and  ||  8  -  0q||  <  p) 

(  I)  ||  being  the  Euclidean  norm  on  Q  )  ,  and  let  Q*(0q,p)  be 
its  image  under  the  transformation  (5.8).  For  0  c  P  let 

(5.13)  n  (0)  =  lim  {  A[ft*  (0 ,  p)  ]  /  A[P(0,p)]} 

p-K) 

where  A  denotes  area.  Then  by  Theorem  VIII  of  Wald  (1943)  the 
likelihood  ratio  test  of  (0.5)  is  asymptotically  optimal  in  the 
sense  that  it 


(a)  has  asymptotically  best  average  power  with  respect  to 
the  weight  function  n(9)  and  the  family  of  surfaces 

(5.14)  $  =  {S(o),c)  :  to  c  fig  ,  c  >  0}  ; 

(b)  has  asymptotically  best  constant  power  on  the  surfaces 

in  S  ; 

and 

(c)  is  an  asymptotically  most  stringent  test. 

By  Corollary  4.3.,  with  the  score-generating  function 
1 1>  =  <p£  the  proposed  rank-order  test  is  asymptotically  power- 
equivalent  to  the  Wald-optimal  likelihood  ratio  test.  Thus  if 
the  underlying  distribution  F  is  logistic,  then  the  QN~test 
using  Wilcoxon-type  scores  is  asymptotically  optimal;  and  if 
F  is  normal,  then  the  normal-scores  rank-order  test  is  asymp¬ 
totically  optimal. 
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